Systematic enumeration of crystalline networks with only sp^ configuration in cubic 

lattices 
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Systematic enumeration of crystalline networks with some special topological characters is of con- 
siderable interest in both mathematics and crystallography. Based on the restriction of lattice in 
cubic and inequivalent nodes not exceeding three, a simple method is proposed for systematic search- 
ing for three-dimensional crystalline networks with only sp^-configuration nodes (C-sp'^-TDTCNs). 
We systematically scan the cubic space groups from No. 195 to No. 230 and find many C-sp^-TDTCNs 
including all the previously proposed cubic ones. These C-sp^-TDTCNs are topologically intrigu- 
ing and can be considered as good templates for searching carbon crystals with novel properties, 
predicting high pressure phases of element nitrogen and designing three-dimensional hydrocarbon 
crystals. Structure optimizations are considered by regrading these C-sp^-TDTCNs as carbon crys- 
tals and the corresponding energetic stability of these carbon crystals are evaluated, using the the 
density functional theory (DFT) based first-principle calculations. Our results are of wide interests 
in mathematics, condensed physics, crystallography and material science. 

PACS numbers: 61.50.Ks, 61.66.Bi, 62.50. -p, 63.20. D- 
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INTRODUCTION 

Systematic enumeration of crystalline networks with 
some special topological characters is of considerable 
interest in both mathematics and crystallography. With 
particular importance, the three-dimensional (TD) four- 
and/or three-connected networks relevant to a wide 
range of systems (especially to the crystal structures, 
such as zeolitic and the well-known diamond) have been 
systematically investigated [H-fllj. TD three-connected 
networks (TDTCNs) have also been considered[ili,|i,i| 
as topologically intriguing structures in mathematics and 
are also significant in crystallography. The crystalline 
ground state of carbon, namely graphite, contains carbon 
atoms with perfect sp^-configuration (the ground state of 
the binary compound of boron and nitrogen, hexagonal 
boron nitride, also prefers such a sp^-conflguration). But 
the graphite structure is not an exact TD network. Many 
efforts have been paid on predicting carbon crystals with 
sp^-configuration in exact TD networks with the belief 
that graphite is not the only way to fill the TD space 
-2-like 
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There are 



with only trigonal sp^-like atoms 
many TDTCNs d i, S H satisfying such' a requirement 
and some of them have been investigated as carbon 
crystals 12 - 3I] showing novel physical properties. The 
second fact implying the signification of TDTCNs (es- 
pecially of those with sp^-configuration, sp^-TDTCNs) 
is about nitrogen. At ambient condition, element 
nitrogen forms the only gas phase of molecular N2 with 
nitrogen-nitrogen triple bond. With the progress in high 
pressure technology, scientists realized a phase transition 
from molecular nitrogen (N2) to an incompressible 



single-bonded cubic gauche phase (cg-N) [40|, |41| . In such 
a novel solid phase, nitrogen atoms form nontraditional 
sp^-hybridization where three electrons bond to three 
neighbours (fill three sp^-hybridized orbits) forming a 
distorted sp^-configuration and the lone pair electrons 
fill the fourth sp^-hybridized orbit. Many followed 
theoretical predictions (42| - |47| also show that high 
pressure nitrogen prefers the TDTCNs with distorted 
sp^-configurations. Thirdly, sp^-TDTCNs are also 
good templates for designing TD hydrocarbon crystals. 
Many efforts have been paid in this direction(48-53']. 
For example, K4-carbon 3^ 37| to K4-hydrocarbon[53], 
graphene jH^I to graphane|48 and graphite to graphane 
crystals [49j. 

Predictions of the TDTCNs can be found in many 
previous works d i, i, i] and the predictions of non- 
graphite TD sp^-hybridized carbon crystals, according 
to our knowledge, were started in 1983 by Hoffmann 
(bct4-carbon) [l^]. Since 1983, many scientists have 
reported their predictions on TD sp^-hybridized car- 
bon crystals, such as H6-carbon in 1990[13], 6.8^D 
(polybenzene) and 6.8^? in 1992*2^, 6.8^0 and P56 in 
1993 25J, carbon forms in 1995^], C20-carbon in 1997 
;27] and later in 1998j2i, FCC-(C28)2, FCC-(C36)2 and 
FCC-(C4fl]2 in 1998 ^23^^ R6-carbon and BCT8-carbon 
in 1998[23, FCC-(C64l2 in 2003(23, K4-carbon in 
2008(11 and 2009^36, sJL C152 and C200 in 2010(2lj, 
TD tubular carbons [Ml3i,j34 Jpjsl and Schwarzite 
carbons (P, D, G surface) (TMI |25| . By combining 
graph theory with quantum mechanics, Winkler et al. 
systematically predicted (sol. (sij 14 TD sp^-hybridized 
carbon crystals with up to six atoms per unit cell. Their 
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results include previously proposed bct4-carbon 12|, 
H6-carbon[13j, R6-carbon[29[ and K4-carbon |35l - l37 |. 
Very recently, with the idea of substitution mentioned 
by Sheng et al (T-carbon) [55*1 , we have also successfully 
constructed [H^ a simple TD sp^-hybridized carbon crys- 
tal (sp^-diamond) with intriguing configuration, which 
belongs to the same space group of diamond and has only 
one inequivalent atomic position. However, the number 
of the TDTCNs is infinite. We can not enumerate all of 
them without any restrictions. In this paper, we report 
a partial solution to this problem, based on a restriction 
of lattice in cubic and inequivalent atomic positions 
not exceeding three. With such a restriction and the 
group theory, a systematic (not a global) enumeration of 
cubic TDTCNs with sp^-like topological configurations 
(C-sp^-TDTCNs) by designed-computational-program 
becomes possible. We have systematically scanned all 
the cubic space groups from No. 195 to No. 230 and 
found many C-sp^-TDTCNs including all the previously 
proposed cubic ones. 



Choose a cubic space group and set the criterion for identifying structures as 
C-sp^-TDTCN [R (e_^.^, 6^^^^) & (R___^__, R__^^^)] 

Index A: Input random position or positions for inequivalent atom or bond or angle 

Initial structure constructing: Operate the inputted positon or positions to 
create a structure and record all the atomic positions in the cubic lattice 

Rationality evaluation if it is a TDTCN (R^^^) 

Caculate the bond angles and length for each inequivalent atoms and find out the 9^^,^ , 

^niax' ^miti ^""^ ^Lx ^^'^ StrUCtUrC. 

If 0\ <0^^^_^^ or 0|^3^>0^^^^7 give up and go to A to create a new structure 

If R^^<R^^^ and I^^^^R^^^. give up and go to A to create a new structure 

If R^^^<R_^^^^ and R^„<R„,„ 

If R^ ^^*R /R^ > R , give up and go to A to create a new structure 
lfR^^*R /R^ <R , record and enlarge the cell with ratio of R /R^ 

If R^. >R^,^ and R^„>R 

If R\ *R /R\ <R , give up and go to A to create a new structure 
^Lin*'^max^'^lax^'^n.in' ^'^'^^^'d and rcducc thc ccU with ratio of R„,^yR-j„^ 

If R* >R and R' <R , record the structure with the initial cell 



FIG. 1: Procedure for systematic searching for C-sp^- 
TDTCNs 



METHOD 

Before we introduce the method of our systematic enu- 
meration of C-sp^-TDTCNs, four simple but very useful 
rules are discussed. Rule one, if a C-TDTCN contains 
only one inequivalent atomic position, the three neigh- 
bours for every atoms are themselves; Rule two, if a 
C-TDTCN contains only two inequivalent atomic posi- 
tions, its two inequivalent atoms form a chemical bond. 
Rule three, if a C-TDTCN contains only three inequiv- 
alent atomic positions, its three inequivalent atoms form 
an angle. Rule four, variation of lattice constant for any 
cubic crystal does not affect its bond angles. Based on 
these four simple rules, we can systematically search C- 
sp^-TDTCNs with one, two and three inequivalent atoms 
by random putting an atom, a bond and an angle, respec- 
tively, in a given cubic lattice in a special space group 
with a defined sp^-criterion. Such an idea enhances the 
probability of generating C-sp^-TDTCNs in comparison 
with a totally random method. A simple procedure im- 
plementing such an idea is shown in Figl. Predictions of 
C-sp^-TDTCNs containing more than three inequivalent 
atomic positions or belonging to non-cubic lattice become 
more complicated (that can also be solved in principle) 
and are not considered in our present work. 

Based on the procedure shown Figl, one can sys- 
tematically search C-sp^-TDTCNs in the following four 
simple steps. I), Choose a cubic lattice with desig- 
nated lattice constant and space group and set a sp^- 
criterion for identifying structures as C-sp^-TDTCNs, 
namely the cutoff radius (Rcut) for the fourth neigh- 
bour, the minima and maxima of bond lengths (Rmin, 
Umax) and bond angles {9min, Qraax) in C-sp^-TDTCNs. 



II), Randomly input the position or positions for the 
inequivalent atom or bond or angle in the chosen lat- 
tice and operate them by corresponding symmetry op- 
erators to construct a testing-crystal. Ill), Calculate 
the bond lengths and bond angles for each inequiva- 
lent atoms in thc created testing-crystal and find out 
the minima (R^j„, 6*^^) and maxima (R^„^, 6*^^^) for 
bond lengthes and bond angles, respectively. IV), Eval- 
uate the rationality of the testing-crystal according to 
the following five conditions if it can be considered as a 
three-connected network (The requirement of Rc„t is sat- 
isfied): a) 9^rnin<(^rmn Or Of^^^>0max indicates that such 
a structure dose not satisfy the requirements of bond 
angle. We give up it directly because further cell en- 
larging or reducing can not change the angles according 
to the rule four mentioned before, b) R^j„<Rmm and 
^max^^^ax indicate that such a structure dose not sat- 
isfy both the minimum and maximum requirements of 
bond length. We give up it directly because further cell 
enlarging (reducing) can not enlarge (reduce) the R^j„ 
(R-tnax) and reduce (enlarge) the R^„^ (R™„), making 
the requirements of bond length satisfied, c) R^j„<Rmin 
and R^^2.<Rmaa; indicate that such a structure dose not 
satisfy the minimum bond length requirement. We then 
enlarge the crystal cell with ratio of Rmin/^min to make 
the crystal satisfying the minimum bond length require- 
ment. If the structure is still a three-connected one 
and the maximum bond length requirement is satisfied 
(namely, Rf^^^*Rmin/'R"rnin<^rnax), wc rccord the en- 
larged structure (otherwise give it up), d) R^j„>Rmin 
and Rf„(j2,>Rmaa: indicate that such a structure dose not 
satisfy the maximum bond length requirement. We then 
reduce the crystal cell with ratio of Rmax/Rfnax make 
the crystal satisfy the maximum bond length require- 
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224-24h-24k-24k(D) 221-48n-48n-24k(P) 217-48h-48h-48h(G) 204-48h-24g(G) 



FIG. 2: Some selected crystalline views of new C-sp^- 
TDTCNs found in present work. Different colors indicate 
inequivalent atomic positions and D, G, P indicate the sur- 
face types. 



ment. If the reduced structure is still three-connected one 
and the minimum bond length requirement is still satis- 
fied (namely, R,''„i„*R max/^max^^miri) ^ ^6 rCCOrd the 
reduced structure (otherwise give it up), e) R^j„>Rmi„ 
and R^^^<Rmax indicate that such a structure directly 
satisfies both maximum and minimum bond length re- 
quirements and we record it directly. 



RESULTS AND DISCUSSION 

In our present work we systematically enumerate the 
C-sp^-TDTCNs by the method as discussed above. We 
name all the C-sp^-TDTCNs with their space group num- 
ber and their corresponding WyckofF positions. For ex- 
ample, the previously proposed 6.8^d[23|, C20-carbon 
11 Q and P192 [i3(P8[il), containing one, two and 
three inequivalent atoms, are named as 224-241, 225- 
48h-32f and 229-961-48k-48k, respectively. With such a 
nomenclature, we can know the fundamental informa- 
tion of a structure, which are useful for identifying if it 
has been proposed previously and excluding repeating 
structures. Based on such a nomenclature, some C-sp^- 
TDTCNs possess the same name due to their common- 
ness in both space group and Wyckoff positions. We pro- 
vide their structural characters to distinguish them. We 
have systematically scanned all the cubic space groups 
from No. 195 to No. 230 based on the requirements of cut- 
off radius of 2.0 A, bond angles distributing in 109°-135° 
and bond lengthes within the range from 1.3 A to 1.7 A 
and found many structurally intriguing C-sp^-TDTCNs 
including all the previously proposed cubic ones (some 
of them are selected to show in Fig2). Especially, C-sp^- 
TDTCNs can be found in every space group but some of 
them degenerate to the ones with relatively higher sym- 
metry. 

The summarized lattice constants and inequivalent 



atomic positions for these C-sp^-TDTCNs are based on 
the post-processes of structure optimizations by regard- 
ing these C-sp^-TDTCNs as sp^-hybridized carbon crys- 
tals. We have considered further energy evaluation for 
these C-sp^-TDTCNs. In the following, we introduce our 
results of C-sp^-TDTCNs templates for sp^-like carbon 
crystals, high pressure crystals of element nitrogen and 
TD hydrocarbon crystals. 



C-sp -TDTCNs with one inequivalent position 

In our searching process, only seven C-sp^-TDTCNs 
containing one inequivalent atomic position were found. 
Although our searching process is non-global, we still 
believe that there are only seven C-sp^-TDTCNs with 
one inequivalent atomic position under the given require- 
ments of cutoff radius, bond lengthes and bond angles. 
We identify that six of them are the previously proposed 
6.8^0 (No.23O)0,j2l], interpenetrated-K4 (No.230, we 
refer it as IK4)[lli76.82p (No.229) [11], sp^-diamond 
(No.227)l56l, 6.8^0 (No.224)El and K4 (No.214)[3ii7| 
crystals. Their crystalline views (fully optimized car- 
bon crystals with first-principles method with calcula- 
tion details provided in the supplementary information) 
are shown in FigSl. (a), (b), (c), (d), (e) and (f), re- 
spectively, in the supplementary information. We name 
these C-sp2-TDTCNs as 230-96h, 230-16b, 229-48k, 227- 
96g, 224-241 and 214-8a, respectively, according to their 
space group numbers and Wyckoff positions. Their cor- 
responding lattice constants and inequivalent atomic po- 
sitions are summarized in Table SI in the supplementary 
information. The seventh one belonging to 206 space 
group has not been reported as a carbon crystal before 
but was proposed as the 3/6/c5 configuration by E. Koch 



et. al in 1995[6| (The structure of sp'^-diamond[56j was 
also proposed as 3/6/c3 by E. Koch 1995[^). Its per- 
spective views of both crystalline and primitive cells are 
shown in FigSl (g) and (h), respectively. From its crys- 
talline view, we can see that it contains eight intercon- 
nected benzene rings in its cubic cell similar to that of 
229-48k but in different connecting manner. We name it 
as 206-48e according to its space group number and the 
fact that it contains one inequivalent atomic position at 
the Wyckoff position of 48e (0.888, 0.985, 0.155). 

We have also found, in some other space group, many 
C-sp^-TDTCNs with only one inequivalent atomic posi- 
tion, which degenerate to the seven ones with relatively 
higher symmetry. For example, we found four equivalents 
for 224-241 (6.8^0) in space groups 228 (228-192h), 227 
(227-192i), 215 (215-24j) and 201 (201-24h), respectively. 
More equivalent relations among these C-sp^-TDTCNs 
with only one inequivalent atomic position are shown in 
Table SI. Furthermore, these seven inequivalent C-sp^- 
TDTCNs can also be found in the following two con- 
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TABLE I: The nominations, characters or previous names (Ch/Pn), lattice constant (LC:A), inequivalent positions (POS), and 
cohesive energy (Ecoh: eV/atom) of C-sp^-TDTCNs with two inequivalent atomic positions 



C-sp2-TDTCNs Ch/Pn 



LC 



POSl 



POS2 



Ecoh 



230-96h-48g 
229-961-961 
227-192i-192i 
227-96g-32e 
227-192i-96g 
226-192j-64g 
225-48h-32f 
225-48g-32f 
225-1921-96k 
225-1921-96j 
224-481-481 
224-12g-8e 
221-24m-241 
221-48n-241 
221-48n-24k 
220-24d-16c 
217-12e-8c 
215-6g-4e 
215-24j-24j 
214-481-481 
212-24e-8c 
205-24d-8c 
204-48h-24g 
200-12j-8f 
195-6g-4o 



New 

New 
D96-D8 
D32-D8 
D72-D9 

New 
C20f27, 28] 

New 

New 

New 
D96-D8 

New 
P48-P9 
P72-P8 
P72-P8 

New 
Cage-N10[47] 

New 

New 
G48-G8 

New 

New 
G36-G7 

New 

New 



10.341 
12.534 
18.048 
11.203 
15.605 
13.227 
9.146 
9.572 
15.146 
14.789 
9.168 
5.342 
8.197 
9.464 
9.462 
6.808 
4.979 
4.904 
8.723 
10.532 
7.705 
6.391 
9.365 
5.451 
4.1303 



■0.144, 
■0.389, 
0.339 
0.834 
0.130 
0.399 
■0.139 
0.250 
0.596 
■0.283 
0.728 
0.000 
0.083 
0.840 
0.339 
0.234 
0.648 
0.500 
0.076 
0.347, 
0.803 
0.410 
■0.133 
0.241 
0.000 



■0.326 
0.231 
■0.716 
■0.166 
0.431 
■0.447, 
0.139 
0.430 
■0.839 
0.352 
0.837, 
0.000 
0.353 
0.371 
0.234 
■0.500 
0.000 
0.500 
0.458 
■0.640 
0.473 
0.390 
■0.342 
0.369 
0.500 



0.776) 
■0.311) 
■0.606) 
0.506) 
■0.749) 
■0.796) 
0.139) 
■0.250) 
■0.728) 
0.418) 
0.947) 
0.362) 
0.353) 
0.734) 
0.129) 
0.250) 
0.000) 
0.140) 
0.347) 
0.196) 
0.305) 
0.672) 
0.274) 
0.000) 
0.322) 



■0.301, 
0.054, 
0.811 
0.387, 
■0.281, 
0.621 
0.000 
0.850 
0.455 
■0.224 
0.323 
0.358 
0.500 
0.072 
0.000 
■0.412 
0.324, 
0.693 
0.677, 
0.007, 
0.715, 
0.258 
0.285 
0.223 
0.241 



■0.875 
■0.222 
■0.412 
0.613 
0.925 
■0.621 
■0.198 
0.650 
■0.307, 
0.500 
0.572 
0.642 
0.838 
0.500 
0.572 
■0.412 
■0.324, 
0.307, 
0.789 
■0.669 
0.715 
0.242 
0.000 
0.777 
0.759 



0.551) 
■0.629) 
■0.534) 
0.613) 
■0.781) 
■0.379) 
0.198) 
■0.650) 
■0.307) 
0.341) 
0.876) 
0.642) 
0.278) 
0.726) 
0.774) 
0.588) 
■0.324) 
0.307) 
0.561) 
0.492) 
0.715) 
0.758) 
0.572) 
0.777) 
0.759) 



-6.824 
-7.278 
-7.443 

■7.191[5r| 
-7.253 
-7.191 
-6.878 
-6.819 
-7.219 
-7.289 
-7.316 

■6.721[5r| 
-7.124 
-7.525 
-7.521 
-7.118 
-6.249 
-6.365 
-7.188 
-7.039 

■6.972f57| 
-6.967 
-7.517 
-7.067 
-6.366 



ditions of considering about two and three inequivalent 
atomic positions in a cubic cell. See the details of equiva- 
lent relations in all the C-sp^-TDTCNs in Table SII and 
SIII in the supplementary information. 



C-sp^-TDTCNs with two inequivalent positions 

When we consider the situation of C-sp^-TDTCNs 
containing two inequivalent atomic positions, we have 
found 25 inequivalent C-sp^-TDTCNs. The perspective 
crystalline views of these 25 inequivalent C-sp^-TDTCNs 
are shown in FigS2 in the supplementary information 
and their corresponding lattice constants and inequiv- 
alent atomic positions are summarized in Table I. 
Most of the 25 C-sp^-TDTCNs are new except that 
the C20-carbon (225-48h-32f) has been proposed as a 
carbon crystaljH [lland the cage-NlO (217-12e-8c) 
has been proposed [47| as a high pressure nitrogen 
phase. Some of them can be regarded as new Schwarzite 
surfaces (indicated in Table I), such as 227-192i-192i (D 
surface with 96 atomic positions in its primitive cell and 
octagon rings, D96-D8), 221-24m-241 (P surface with 
48 atomic positions in its primitive cell and nonagon 
rings, P48-P9) and 204-48h-24g (G surface with 36 



atomic positions in its primitive cell and heptagon 
rings, G36-G7). These C-sp^-TDTCNs distribute in 
space groups of 230, 229, 227, 226, 225, 224, 221, 220, 
217, 215, 214, 212, 205, 204, 200, and 195. Some of 
them can also be found in other space groups or in 
the following condition of C-sp^-TDTCNs containing 
three atomic positions (see the equivalent relations in 
Table SII and SIII in the supplementary information). 
In addition, we have also found some C-sp^-TDTCNs 
with two inequivalent atomic positions degenerating to 
those with only one inequivalent atomic position. Such 
as the one named as 228-96g-96g is equivalent to 224-241 
(6.8^D). More examples can be found in Table SII in the 
supplementary information. 



C-TDTCNs with three inequivalent positions 

There are more C-sp^-TDTCNs containing three 
inequivalent atomic positions and the increase of the 
inequivalent atomic positions enlarges the state space 
and reduces the probability of meeting the special C-sp^- 
TDTCNs in our searching process. Fortunately, we have 
found 48 of them including some Schwarzite type surface 
structures (indicated in Table II), such as the G type 



5 



TABLE II: Nominations, characters or previous names (Ch/Pn), Lattice constant (LC:A), Inequivalent positions (POS), Co- 
hesive energy (Ecoh: eV/atom) of C-sp^-TDTCNs with three inequivalent atomic positions 



C-sp^-TDTCNs 


Ch/Pn 


LC 




POSl 








POS2 










POS2 






Ecoh 


9*^0 QfiVi QfiVi QfiVi 


ni 44 n8 


15 944 


-0.054 


-0.491 


0.569) 


(-0 


152 


-0 


559, 





661) 


(-0 


099 


-0 


488, 





644) 


-7 


531 


2'^n-4Sp--4Sf-1 fih 


New 


9.882 


-0.021 


-0.375 


0.271) 


( 


000 


-0 


825, 





451) 


(-0 


125 


-0 


375, 





375) 


-7 


034 


QfiVi QfiVi 4Sf 




12.038 


-0.035 


-0.342 


0.501) 


(-0 


134 


-0 


400, 





531) 


( 


000 


-0 


250, 





564) 


-7 


283 


22Q-Qfil-4Sk-1 fif 

^^tj iy\Jl. rtijl\. J-\JL 


Nol 


11.844 


0.229 


0.638 


0.059) 


( 


413 





263, 


-0 


263) 


( 


317 





317, 





317) 


-7 


125 


99Q dHi dHV 


New 


11 739 


0.241 


-0.156 


0.000) 


( 


327 


-0 


187, 


-0 


187) 


( 


281 


-0 


111, 


-0 


111) 


-7 


013 


99Q QRl 4Slr zL^^V 

iyUi rlOiV rtO-Ti 


PI Q9 PJ^ri 7 1 Ql 


14.901 


0.278 


-0.047 


0.396) 


( 


289 


-0 


173, 





289) 


f 


324 


-0 


091, 





324) 


-7 


627 


99^^ 1Q91i Qficr '^9h 


New 


14.433 


-0.326 


0.857 


0.542) 


(-0 


375 





698, 





553) 


(-0 


375 





625, 





625) 


-7 


175 


997 QficT Qficr ^9p 


New 


12.652 


0.169 


-0.533 


0.830) 


( 


068 


-0 


288, 





788) 


( 


171 


-0 


171, 





671) 


-7 


207 






17.491 


-0.105 


0.605 


0.801) 


(-0 


244 





858, 





858) 


(-0 


123 





778, 





778) 


-7 


125 


227-1 Q2i-1 Q2i-Qfip- 


7)1 90-7)7 


18.994 


-0.402 


-0.711 


-0.005) 


(-0 


910 


-0 


804, 





357) 


(-0 


725 


-0 


275, 


-0 


030) 


-7 


518 


227-1 Q2i-Qfi!r-'^2p 




15.881 


-0.083 


-0.021 


0.703) 


( 


063 


-0 


326, 





937) 


( 


108 


-0 


392, 





892) 


-7 


324 


227-1 Q2i-1 q2i-1 Q2i 


D144-D8 


23 623 


-0.315 
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230-96h-96h-96h (G144-G8) and P type 229-961-48k-48k 
(P192-P8, proposed before[l3, El). Most of them are 
never mentioned in literatures. Their corresponding 
lattice constants and inequivalent positions are summa- 



rized in Table II and their corresponding perspective 
crystalline views can be found in the supplementary Fig 
S3 and Fig S4. Some of them are identified to be the 
previously proposed FCC-(C36)2l23i] (227-96g-96g-96g), 
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FCC-(C4o)2l23 (227-192i-96§^2eK FCC-(C28)2l2 



(203-96g-96g-32e) and P56-P7|16l (195-12j-12j-6g). 
We notice that many new surface structures, such 
as 227-192i-192i-96g (D120-D7), 227-192i-192i-192i 
(D144-D8), 224-481-481-241 (D144-D8), 221-24m-241-8g 
(P56-P8), 221-48n-48n-24k (P120-P8), 217-48h-48h-48h 
(G72-G8), and 201-24h-24h-8e (D56-D7) have not been 
proposed previously. Some new structures belonging to 
the same space and with the same atomic positions have 
the same name (see 217-48h-48h-48h, 215-24j-24j-24j 
and 206-48e-48e-48e) and we provide additional signs 
reflecting their structural characters for distinguishing 
them. For example, 6*A7^ means the structure with 
characteristic hexagons connecting six heptagons at its 
six armchair edges and 6*Z6^ indicates characteristic 
hexagons connecting six hexagons at its six zigzag 
vertexes. 

From table II, we can see that these C-sp^-TDTCNs 
distribute in space groups of 230, 229, 228, 227, 226, 
225, 224, 223, 221, 217, 217, 216, 215, 214, 212, 207, 
206, 204, 203, 201, 200, 197 and 195. This dose not 
means that C-sp^-TDTCNs can not be found in other 
space groups. For example, we found 227-96g-96g-32e 
in No.210 as 210-96h-96h-32e with relatively lower 
symmetry. In addition, we have also found some 
C-sp^-TDTCNs with three inequivalent atomic positions 
degenerating to those with one or two inequivalent 
atomic positions. Such as 213-8c-4a-4b to 230-16b (IK4) 
and 220-48e-48e-48e to 230-96h-48g. More examples can 
be found in Table SIII in the supplementary information. 



Additional discussions 

All the post-processes of structure optimizations 
for these C-sp^-TDTCNs are performed with density 
functional theory (DFT) based first-principles method 
by regarding these C-sp^-TDTCNs as sp^-hybridized 
carbon crystals with calculation details provided in the 
supplementary information. We should point out that, 
after the optimization, a few C-sp^-TDTCNs do not sat- 
isfy the requirements of C-sp^-TDTCNs in bond angles, 
bond lengths or neighbours, which are understandable 
and acceptable. We have also considered energy evalua- 
tion of these sp^-hybridized carbon crystals but have not 
performed systematic investigation on the dynamical 
stability due to the large computational cost. Our results 
show that some of them are energetically viable, such 
as 230-96h-96h-96h (G144-G8), 229-961-48k-48k (P192- 
P8), 227-192i-192i-96g (D120-D7), 224-24h-24k-24k 
(D72-D8), 221-48n-48n-24k (P120-P8), 217-48h-48h-48h 
(G72-G8), 221-48n-24k (P72-P8), 221-48n-241 (P72-P8), 
204-48h-24g (G36-G7), and 224-241 (6.8^0) which pos- 
sess remarkable energetic stability comparable to that of 
diamond (-7.623 eV/atom). Most of them are more sta- 



ble than the experimentally produced graphdiyne (-7.06 
eV/atom). In addition, we want to point out that some 
of them rr iay be dynamically unstable like K4-carbon 
(214-16b) [37[, and some of them such as sp^-diamond 
and 6.8^D have been confirmed dynamically stablejs^. 



CONCLUSION 

Based on a simple restriction of lattice in cubic and 
inequivalent atomic positions not exceeding three, we 
have introduced a simple method for systematic search- 
ing of C-sp^-TDTCNs in all the cubic space groups from 
No. 195 to No. 230 and found many topological intriguing 
C-sp2-TDTCNs. These C-sp^-TDTCNs are named 
according to their space group numbers and Wyckoff po- 
sitions, which can provide fundamental crystallographic 
information for excluding future repeating structures. 
All of these C-sp^-TDTCNs can be considered as good 
templates for searching for carbon crystals with novel 
properties, predicting high pressure phases of element 
nitrogen and designing TD hydrocarbon crystals. We 
believe that our results are of wide interests in both 
mathematics and crystallography. 
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We optimized these C-sp^-TDTCNs through regarding them as carbon crystals. Both the lattice 
constant and atomic positions of these C-sp^-TDTCNs were fully optimized until the residual forces 
on each carbon atom to be less than 0.02 eV/A. All the calculations of structure optimizations 
and energetic stability evaluations were performed using the density functional theory based VASP 
code [l|] with the projected augmented wave (PAW) potential Q]. The exchange and correlation are 
approximated by general gradient approximation (GGA) developed by Perdew et al. [^. The wave 
functions for all systems are expanded by plane-wave functions with cutoff energy of 500 eV. The 
Brillouin zone sample meshes based on the Monkhorst-Pack scheme are set to be denser enough 
(less than 0.25 1/A) to ensure the accuracy of our calculations. 

PACS numbers; 61.50.Ks, 61.66.Bi, 62.50. -p, 63.20. D- 
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TABLE I: The nominations, relations, lattice constant (LC:A), inequivalent positions (POS), and cohesive energy (Ecoh: 
eV/atom) of C-sp^-TDTCNs with only one inequivalent atomic position 



C-TDTCNs 


R,clcitiion 


LC 




POS 




Ecoh 


230- 96h 


6.8^G 


9.587 


( 0.117, 


0.548, 


0.424) 


-7.391 


230-1 fib 


J 4- car b o n 


4.426 


(-0.125 


-0.375 


0.375) 


-6.346 


229-48k 


6.8^P 


7.843 


( 0.689 


0.689 


0.913) 


-7.364 


227-96 e- 


sp^-diSiinoiici 


9.667 


(-0.049 


0.549 


-0.726) 


-7.179 


224-241 


6.8^D 


6.095 


( 0.250 


0.087 


0.589) 


-7.585 


214-8a 


K4-c3.rb)ori 


4.126 


(-0.125 


-0.375 


0.375) 


-6.529 




3 /fi /r5 


7.505 


( 0.345 


-0.111 


0.015) 


-6.798 




6.8^D 


12.191 


( 0.207 


0.043 


0.875) 


-7.585 


997-1 Q9i 


6.8^D 


12.191 


( 0.625 


0.793 


0.543) 


-7.585 


99'^-A«l 

Zi^O rtOl 


6.8^P 


7.843 


( 0.311 


0.087 


0.689) 


-7.364 


999 4Si 


g g2p 


7.843 


( 0.311 


0.087 


0.689) 


-7.364 


990 1 fir 




4.426 


(-0.125 


-0.375 


0.375) 


-6.346 


215-24j 


6.8^D 


6.095 


( 0.250 


0.087 


0.589) 


-7.585 


2 13- 8c 


Tv 4— r* a r\ n 


4.126 


(-0.125 


-0.375 


0.375) 


-6.529 


212-8C 


K4- carbon 


4.126 


( 0.375, 


0.375, 


0.375) 


-6.529 


211-48j 


6.82p 


7.843 


( 0.689, 


0.689, 


0.913) 


-7.364 


210-96h 


sp^-diamond 


9.667 


(-0.049, 


0.549, 


-0.726) 


-7.179 


206-16C 


IK4- carbon 


4.426 


(-0.125, 


-0.375, 


0.375) 


-6.346 


204-48h 


6.82p 


7.843 


( 0.689, 


0.689, 


0.913) 


-7.364 


203-96g 


sp^-diamond 


9.667 


(-0.049, 


0.549, 


-0.726) 


-7.179 


201-24h 


6.82D 


6.095 


( 0.250, 


0.087, 


0.589) 


-7.585 


199-8a 


K4-caj:bon 


4.126 


(-0.125, 


-0.375, 


0.375) 


-6.529 
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TABLE II: The nominations, relations, lattice constant (LC:A), inequivalent positions (POS), and cohesive energy (Ecoh: 
eV/atom) of C-sp^-TDTCNs with two inequivalent atomic positions 



C-TDTCNs 


R clcitioii 


LC 






POSl 










POS2 






Ecoh 


22S-1 Q2h-1 Q2h 


g 82p 


15.685 


( 


043 





155, 





845) 


f 


043 


0.155 





655) 


-7 


364 


yug yug 




12.191 


( 


875 





289, 


Q 


539) 


f 


291 


0.375 


Q 


041) 


-7 


585 


997-QRVi QfiVi 


6.8^D 


12.191 


f 


467 


-0 


375, 


Q 


783) 


f 


462 


0.712. 


Q 


375) 


-7 


585 


99fi-1Q9i 1Q9i 


6.8^P 


15.685 


( 


155 





957, 


Q 


845) 


f 


207 


0.905 


Q 


905) 


-7 


364 


991 94m 94m 




7.843 


f 


311 





087, 


Q 


689) 


f 


413 


0.189 


Q 


811) 


-7 


364 


Zi^Lf T:OtJ ^OC 


g gZQ 


9.587 


f-0 


798 


-0 


363, 


-0 


674) 


f-0 


617 


-0.547 


-0 


076) 


-7 


391 


91 Q-QRh-QfiVi 


6.8^D 


12.191 


( 


125 





043, 


Q 


793) 


f 


207 


0.043 


Q 


875) 


-7 


585 


91 Q 4Scr ^9o 


21 ^ f^o- 4o 


9.808 


( 


250 





070, 


Q 


750) 


f 


154 


0.154 


Q 


846) 


-6 


365 


91 M-94i-94i 


6.8^P 


7.843 


( 


087 


-0 


311, 





311) 


f-0 


087 


-0.311 





311) 


-7 


364 


218-12f-8p 


21 7-1 2p-8r 


4.979 


( 


648 





000, 





000) 


f 


324 


-0.324 


-0 


324) 


-6 


249 


21 7-24fr-24(>- 


g,g2p 


7.843 


( 


087 


-0 


311, 





311) 


f-0 


087 


-0.311 





311) 


-7 


364 


21 fi-4Sh-4Sh 


— rl 1 p m nn H 

IXlCtlllL'llLl 


9.667 


(-0 


049 





549, 


-0 


726) 


f-0 


201 


0.476 


-0 


299) 


-7 


179 


91 fi-Qfii-QRi 


6.8^D 


12.191 


( 


043 





125, 


Q 


293) 


f 


043 


0.125 





793) 


-7 


585 




T rv 4_i^Q r^n 


4.426 


f-0 


125 


-0 


375, 


Q 


375) 


f 


125 


-0.625. 


Q 


625) 


-6 


346 


914-48i-48i 

tOl t:01 


6.8^G 


9.598 


f 


452 





424, 


Q 


617) 


f 


576 


0.383 


Q 


548) 


-7 


391 


91 4a A\^ 


T\ /I /^ovr^i^Vl 
XVi L-Ctl UUii 


4.126 


( 


625 





625, 


Q 


625) 


f 


125 


-0.625 


Q 


625) 


-6 


529 


91 9-4a-4l-> 


xvt L-di uuii 


4.126 


f 


875 





625, 


Q 


375) 


f 


625 


0.625 


Q 


625) 


-6 


529 


21 n-Qfih-S2p 


227-Qfih-S2p 


11.203 


I 


834 


-0 


166, 





506) 


f 


387 


0.613 





613) 


-7 


191 


2nQ-4Sh-S2f 

^\JtJ rrOll Oirfi- 


C20 


9.146 


(-0 


139 


-0 


139, 





1 '^91 


f 


000, 


-0.198 


.0 


198) 


-6 


878 


209-48i-32f 


225-48g-32f 


9.572 


( 


250 





430, 


-0 


250) 


f 


850 


0.650 


-0 


650) 


-6 


819 


2n8-24m-24Tn 


6.8^P 


7.843 


( 


087 


-0 


311, 





311) 


f 


587 


-0.189 





189) 


-7 


364 


2(18-1 21-1 21 

X Zl X ^1 


6.8^D 


6.095 


( 


250, 





913, 


Q 


413) 


f 


250, 


0.587, 


Q 


087) 


-7 


585 


9n7-94V 94V 


g_g2p 


7.843 


( 


087 


-0 


311, 


Q 


311) 


f 


587 


-0.189 


Q 


189) 


-7 


364 


9n^-><r'->^r 
^ u u oc oc 


TT^4_f aT"V\r»n 


4.426 


f-0 


125 


-0 


375, 





375) 


f-0 


125 


-0.625 





125) 


-6 


346 


9nV94H-94fl 




7.505 


f 


595 





235, 


-0 


361) 


f 


095 


0.735 





139) 


-6 


798 


9rm Qfio- ^9p 

zuo yug ozrc 


297 Qficr '^9p 


11.203 


f 


834 


-0 


166, 


Q 


506) 


f 


387 


0.613 


Q 


613) 


-7 


191 


909 4Kli '^9f 

ZiUii rtOli Oiii 


C20 


9.146 


f-0 


139 


n 
-u 


1 "^Q 

-L o y , 


Q 


139) 


f 


000 


-0.198 


Q 


198) 


-o 


878 


202-48g-32f 


225-48g-32f 


9.572 


( 


250, 





430, 


-0 


250) 


f 


850, 


0.650, 


-0 


650) 


-6 


819 


201-24h-24h 


6.82p 


7.843 


f 


087, 


-0 


311, 





311) 


f 


587, 


-0.189, 





189) 


-7 


364 


201-12f-8c 


224-12g-8e 


5.342 


f 


000, 





000, 





362) 


f 


358, 


0.642, 





642) 


-6 


721 


198-4a-4b 


K4- carbon 


4.126 


( 


875, 





625, 





375) 


( 


625, 


0.625, 





625) 


-6 


529 


197-24g-24g 


6.82p 


7.843 


f 


087, 


-0 


311, 





311) 


f-0 


087, 


-0.311, 





311) 


-7 


364 


197-12c-8c 


CagclO-N 


4.979 


( 


648, 





000, 





000) 


f 


324, 


-0.324, 


-0 


324) 


-6 


249 


196-48h-48h 


sp'^-diamond 


9.667 


f-0 


049, 





549, 


-0 


726) 


f-0 


201, 


0.476, 


-0 


299) 


-7 


179 


195-12j-12j 




6.095 


( 


250, 





913, 





413) 


( 


250, 


0.587, 





087) 


-7 


585 
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TABLE III: Nominations, relations, Lattice constant (LC:A), Inequivalent positions (POS), Cohesive energy (Ecoh: eV/atom) 
of C-sp^-TDTCNs with only one inequivalent atomic position 



C-TDTCNs 


Relation 


LC 


POSl 


POS2 


P0S3 


Ecoh 



228-192h-192h-192h 


224-24h-24k-24k 


20.233 


(-0 


343, 





845, 





577) 


( 


440, 





784, 





940) 


( 


501, 





750, 





910) 


-7 


545 


226-192j-192j-192j 


221-48n-241 


18.927 


( 


420, 





315, 





867) 


( 


464, 





250, 





863) 


{ 


367, 





314, 





920) 


-7 


525 


220-48e-48e-48e 


230-96h-48g 


10.340 


(-0 


144, 


-0 


326, 





776) 


(-0 


301, 


-0 


875, 





551) 


{ 


394, 


-0 


526, 





576) 


-6 


824 


216-48h-16e-16e 


C20 


9.125 


( 


547, 





453, 





272) 


( 


373, 





127, 





627) 


( 


610, 





390, 





390) 


-6 


878 


213-8c-4a-4b 


IK4-carbon 


4.426 


(-0 


125, 


-0 


375, 





375) 


( 


375, 





375, 





375) 


( 


125, 


-0 


625, 





625) 


-6 


346 


212-8c-4a-4b 


IK4-carbon 


4.426 


( 


125, 


-0 


625, 





625) 


(-0 


125, 


-0 


375, 





375) 


(-0 


125, 


-0 


625, 





125) 


-6 


346 


210-96h-96h-32e 


227-96g-96g-32e 


12.652 


( 


169, 


-0 


533, 





830) 


{ 


068, 


-0 


288, 





788) 


( 


171, 


-0 


171, 





671) 


-7 


207 


210-96h-96h-96h 


227-96g-96g-96g 


17.491 


(-0 


105, 





605, 





801) 


(-0 


244, 





858, 





858) 


(-0 


123, 





778, 





778) 


-7 


125 


2u9-9dj-9dj-9dj 


225-1921-9Dk 


15.146 


( 


596, 


-0 


839, 


-0 


728) 


( 


455, 


-0 


307, 


-0 


307) 


( 


404, 


-0 


161, 


-0 


272) 


-7 


219 


209-9dj-9dj-32i 


225-9Dk-9Dk-32t 


13.221 


( 


319, 


-0 


319, 





449) 


{ 


377, 


-0 


623, 





744) 


( 


337, 


-0 


163, 





337) 


-7 


076 


2U7-24k-24k-24k 


221-48n-24k 


9.462 


( 


340, 





234, 





129) 


( 


000, 





572, 





774) 


( 


660, 





766, 





871) 


-7 


521 


205-24k-24k-24k 


221-48n-24k 


9.462 


( 


340, 





234, 





129) 


( 


000, 





572, 





774) 


( 


660, 





766, 





871) 


-7 


521 


203-96g-96g-32e 


227-96g-96g-32e 


12.652 


( 


169, 





967, 





330) 


{ 


213, 





713, 





932) 


( 


171, 





829, 





671) 


-7 


207 


203-96g-96g-96g 


227-96g-96g-96g 


17.491 


(-0 


105, 





605, 





801) 


(-0 


244, 





858, 





858) 


(-0 


123, 





778, 





778) 


-7 


125 


202-96j-96j-96j 


225-1921-96k 


15.146 


( 


596, 


-0 


839, 


-0 


728) 


( 


455, 


-0 


307, 


-0 


307) 


( 


339, 


-0 


096, 


-0 


272) 


-7 


219 


202-96i-96i-32f 


225-96k-96k-32f 


13.221 


( 


319, 


-0 


319, 





449) 


( 


377, 


-0 


623, 





744) 


( 


337, 


-0 


163, 





337) 


-7 


076 


202-96i-32f-32f 


200-12j-8f 


10.901 


( 


250, 





864, 





932) 


( 


364, 





864, 





864) 


( 


135, 





864, 





864) 


-7 


067 


201-24h-24h-24h 


204-48h-24g 


9.365 


{ 


500, 





928, 





785) 


{ 


726, 





867, 





658) 


( 


367, 





842, 





774) 


-7 


517 


197-24f-24f-24f 


204-48h-24g 


9.365 


( 


500, 





928, 





785) 


( 


367, 





842, 





774) 


( 


342, 





726, 





867) 


-7 


517 


195-12j-4e-4e 


200-12j-8f 


5.412 


( 


499, 





271, 


-0 


891) 


( 


718, 





718, 


-0 


282) 


( 


285, 


-0 


285, 


-0 


715) 


-7 


067 




FIG. 1: Crystalline views of the seven inequivalent C-sp^-TDTCNs. (a), (b), (c), (d), (e) and (f) are 230-96h (6.8^G), 230-16b 
(IK4), 229-48k (6.8^P), 227-96g (sp^-diamond), 224-241 (6.8^D) and 214-8a (K4), respectively, (g) and (h) are the crystalline 
views of 206-48e (3/6/c5) in crystal cell and primitive cell, respectively. 
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230-96h-96h-96h 230-48g-48f-16b 230-96h-96h-48f 229-961-48k-16f 229-48j-48k-48k 229-961-48k-48k 




228-192h-96g-32b 227-96g-96g-32e 227-96g-96g-96g 227-192i-192i-96g 227-192i-96g-32e 227-192i-192i-192i 




226-96i-64g-48e 225-96k-96k-32f 225-1921-1921-1921 225-1921-1921-96j 224-481-481-241 224-24k-24k-24k 




224-24h-24k-24k 224-481-481-24j 223-481-24k-12h 221-24m-241-8g 221-48n-48n-24k 221-48n-48n-241 



FIG. 3: Crystalline views of the 24 inequivalent C-sp^-TDTCNs (No. 230-221) with three inequivalent positions, where different 
atomic positions are shown in different colors. 
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217-48h-48h-48h-S 217-481i-481i-48h 217-48h-481i-24g 216-96e-48h-16e 216-96e-96e-48h 215-24i-12i-4e 




215-24j-24j-24j 215-24j-24j-24j 214-48i-48i-48i 212-12d-8c-4a 2i2-24e-24e-8c 207-24k-24k-8g 




206-48e-48e-48e 206-48e-48e-48e 206-48e-48e-48e 204-481i-48h-24g 203-96g-96g-32e 201-24h-24h-8e 




200-12k-8i-6h 200-241-12j-8i 200-12k-8i-6f 200-12j-8i-6g 197-24f-24f-8c 195-12j-12j-6g 



FIG. 4: Crystalline views of the 24 inequivalent C-sp^-TDTCNs (No. 216-195) with three inequivalent positions, where different 
atomic positions are shown in different colors. 



